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Abstract 

We construct electrically charged AdSs black hole solutions whose charge, mass and 
boost-parameters vary slowly with the space-time coordinates. From the perspective of 
the dual theory, these are equivalent to hydrodynamic configurations with varying chemical 
potential, temperature and velocity fields. We compute the boundary theory transport 
coefficients associated with a derivative expansion of the energy momentum tensor and R- 
charge current up to second order. In particular, for the current we find a first order transport 
coefficient associated with the vorticity of the fluid. 
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1 Introduction 

The AdS/CFT correspondence has proved to be a useful tool in understanding various 
aspects of strongly coupled gauge theories. In this work we focus on some developments 
which allow one to relate the hydrodynamic regime of the gauge theory to black hole solutions 
in asymptotically AdSs backgrounds. The first work in this direction was carried out in [1] 
where the ratio of the shear viscosity t] to the entropy density s of the M = 4 SU(N) 
super symmetric Yang-Mills theory was computed via the Kubo formula. It was found that 
at strong t' Hooft coupling and in the large N limit 

- = T - - (1) 
s An w 

This value seems to be universal, and applies to a large class of theories which have a 
holographic dual [21 El HI El El [7], see [8] for a review. Finite N corrections to (CD) were 
considered in [91 [TQl El HH UM [131 [7J . The fact that the black hole background allows to 
compute hydrodynamic transport coefficients might be an indication that black holes capture 
the full hydrodynamic behavior of the boundary theory. In [T4j an important step in this 
direction was made: it was shown how to map a hydrodynamic expansion of the boundary 
theory to a gradient expansion in the bulk. In principle, this technique allows to compute all 
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the transport coefficients of conformal fluid dynamics (for gauge theories with an AdS dual). 
Those coefficients which are accessible via linear response theory can also be computed using 
the Kubo formula. For instance, in [151 HSj some of the second order transport coefficients 
were computed this way. 

The method of [2] has been shown to be rather robust and can be applied to black 
holes in various dimensions [TTJ [T8| , and to situations where the black hole metric couples to 
external fields such as the dilaton [19] , implying forced fluid dynamics in the boundary theory. 
The method involves extending known asymptotically AdS^ black hole solutions by allowing 
various parameters of the solution to vary with the space-time coordinates. In this work we 
consider charged AdSs black holes. We show that the known Reissner-Nordstrom charged 
black hole solutions can be extended so that their charge, mass and certain boost parameters 
are slowly varying in the coordinates transverse to the AdS radial direction (henceforth, the 
transverse coordinates). In the dual picture, this corresponds to a hydrodynamic limit of 
the theory where the charged current, energy density, and velocity fields are slowly varying. 

In detail, the bulk theory we have in mind is Einstein-Maxwell gravity with a negative 
cosmological constant and a Chern-Simons term. This is a consistent truncation of IIB 
supergravity on AdS 5 x S 5 and is dual to the strongly coupled, planar limit of the M = 4 
SU(N) supersymmetric- Yang-Mills theory on R 3,1 with a non vanishing chemical potential 
[201 [21]. We will sometimes call this SYM or M = 4 theory for short. More precisely, it is dual 
to a subsector of the M = 4 theory in which a single conserved U(l) current is excited. This 
is the Noether current associated with the diagonal U(l) of the maximal Abelian subgroup 
of the 50(6) i?-symmetry group and it is dual to the bulk U(l) gauge field. More details on 
this truncation can be found in [2~T] . 

Previous computations of the thermodynamic properties of the SYM fluid with a finite 
chemical potential can be found in [211 EUJ [22] where the energy density and equation of 
state have been analyzed. The shear viscosity of this fluid has been computed via the Kubo 
formula in [51 [231 121], in [5] its heat conductivity was analyzed and in [16] some of the 
dispersion relations were computed to second order in a small momentum expansion. In [25] 
one can find a related analysis dealing with M2-branes. In this work, we extend these results 
and compute all second order, linear and non-linear, transport coefficients. We also find a 
first order contribution to the i?-charge current which was not considered in the literature 
so far. 

The rest of this paper is organized as follows: in the subsequent section we discuss 
conformal fluid hydrodynamics, set the notation for the rest of this paper and summarize 
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our field theory results. In section [3] we review the Reissner-Nordstrom AdSs black hole 
solution and rederive the thermodynamic properties of the associated boundary theory. Our 
main computation, extending the Reissner-Nordstrom black hole solution to one with a 
slowly varying charge, mass and boost parameters is done in section HI There, we also 
explain how to compute the associated transport coefficients. We end with a discussion in 
section Some of the details of the analysis are relegated to an appendix. 

Towards the end of this work we learned about [26] which has some overlap with the 
material presented here. 

2 Conformal fluid dynamics (and summary) 

Consider the energy momentum tensor (T^) and a conserved current (J M ) of a conformal 
theory in R 3,1 . In the absence of anomalies one has 

9 M (T^)=0 9 M (J") = 0, (2) 

and the energy momentum tensor is traceless 

(T£) = 0. (3) 

In the hydrodynamic approximation, where the mean free path of the theory £ m f p is smaller 
than the typical inverse momentum scale, both (T^) and (J M ) can be expressed in terms of 
hydrodynamic fields. These are given by the energy density e, the charge density p and the 
velocity field m m (which is normalized such that u^u^ = —1). In the Landau frame, which we 
use in the rest of this work, all the hydrodynamic fields are defined relative to the rest frame 
of a fluid element. The energy density is given by the time-time component of the energy 
momentum tensor in the rest frame of the fluid element, the charge density is given by the 
zero component, J°, of the current in the rest frame of a fluid element, and the velocity of 
a fluid element is defined by the boost parameter needed to bring that fluid element to its 
rest frame. In general, one can exchange the energy density and charge density with the 
temperature T and chemical potential /i, though their explicit functional relation will depend 
on the specific details of the theory. For the M = 4 theory these relations were computed in 
[2~T | |2"2" | 15] and are given in fl2Ual) and (120b|) together with (I2ip . In, for instance, [27J one can 
find such relations when flavored matter is introduced. 

Working in the hydrodynamic regime, where the momentum scale is smaller than the 
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inverse mean free path, implies that the velocity field, energy density and charge density 
vary slowly with the space-time coordinates, i.e. their derivatives are small. For example, 
one has 

\de\ « e/£ mfp . (4) 

In this case, we can expand the energy momentum tensor and current in gradients of the 
hydrodynamic variables. At zero order in such a gradient expansion (meaning a fluid with 
constant energy density, charge density and moving at a fixed velocity), the only current 
which can be constructed from u^, e and p will be proportional to u^. The only symmetric 
traceless tensor one can construct must be proportional to Tj^ + UpU,,. Thus, to leading order 
in gradients, 

(T^u) = ~ (4W M W^ + 7^) , (Jy) = PU V . (5) 

We denote higher order gradient corrections to the energy momentum tensor and current by 



11^ and T M , 



(Tpv) = | (4uj.it,, + + 11^ , (J J) = pu v + T„. (6) 



Working in the Landau frame, we find that u v T v = and uTI^ = 0. 

Following [15] it is possible to construct the form of the corrections T M and 11^,, order 
by order in a derivative expansion. Let us start by evaluating all possible contributions to 
Y„ at first order in a gradient expansion. Overall, there are four possible vectors that can 
be constructed from e, p and u v which are orthogonal to the velocity field and have one 
derivative. These are 

Vl = P M »d v e Vl = P, v d vP Vl = P M v u a d a u v 

Vl = £, , (7) 

where we have defined 

^ = ^ paT u p d a u T , (8) 
and P M „ projects onto the space orthogonal to the velocity field, 

P^ = u ll u v + r} IJU ,. (9) 

reduces to the curl of the velocity in the local rest frame. In the following, we will always 
adorn vectors or tensors involving with a tilde. 



4 



Since energy conservation ([2]) implies that 



Vl = -4eV u 3 , (10) 

we can construct the leading derivative terms in T„ by using combinations of only Vj-, V 2 and 
Vj}. Further, recall that in a conformal theory, a conserved current should transform homo- 
geneously under Weyl transformations (we use the conventions of [15] where r]^ — > e -2 ^?]^). 

is Weyl invariant while and V 2 transform inhomogeneously with the inhomogeneous 
terms given by 

gyi = e ^ P ^4ed u u , bVl = P/3pd u uj. (11) 
However, the linear combination 

Vfty ( 12 ) 

does transform homogeneously. One could have guessed this by noting that the weight of 
e under Weyl rescalings is 4 and the weight of p under Weyl rescalings is 3, so the only 
Weyl invariant combination is e 3 / p 4 . Switching from energy density and charge density to 
temperature T and chemical potential /i, we find that the most general form of T u at first 
order in a derivative expansion is 

t v = -KP v a d a ^ + m v + o(d 2 ) , (13) 

where k = n(p, T) and Q = Q(fi, T) are undetermined first order transport coefficients whose 
explicit form depends on the theory. Our value for k and Q can be found in (|22bl) and (I22cl) 

In principle, the transport coefficient k in ffl3|) can also be calculated via linear response 
theory. In [5] it has been calculated for a different sector of the M = 4 theory in which a 
different U(l) subgroup of the i?-charge current is excited in addition to some of the scalar 
fields. The coefficient Q in (1131) has been inaccessible so far E| 

In [28] a different type of argument has been used to construct the first order terms in 
T^: an entropy current was constructed by hand to have a positive semi-definite divergence 
and from it, the first order corrections to T u were inferred. The Yj, constructed in [2B] differs 
from the one in ffTB"]) by the i v term. It would be interesting to find a corrected form of the 
entropy current which allows for the l v term, perhaps along the lines of |19j . 

Second order contributions to T u can be derived using the same arguments as those 

1 From the bulk point of view, it is the Chern-Simons term in the action (cf. (|26p below) that is responsible 
for having f2 ^ 0. 
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leading to (|T3|) . One may construct all possible Weyl-covariant vectors composed of two 
derivatives which are orthogonal to the velocity field. See for example [2§J [3U1 ED E21 E2J EDS 
[15] for a more elaborate discussion on the construction of Weyl invariant quantities. In this 
work five such terms will be relevant^] 

s? } = ^4 , 4 2 ) = ^a a | , s?) = p^ q (a>- 3 ) , s?) = P^ Q t^- 1 ) 

S« = , (14) 

where we have defined 

<V = Zd^u u) , uj^ = )^P^P° {d x u a - d a u x ) (15) 

and angular brackets denote a traceless projection onto the space orthogonal to u M so that 

A {H = P^ 1 - (A Xa + A aX ) - J— PfMU P^A Xa (16) 

satisfies rf v A^ v ) = and u^A^ = 0. For the i?-charge current of the theory at hand, the 
only non vanishing second order transport coefficients are those associated with the terms 
in ffl4"|) . Thus, up to second order in a derivative expansion, we may write 

T„ = -KP?d a ± + m v + £iS« + 6H?> + CsH® + e 4 sW + |iH«. (17) 

The decomposition of into Weyl invariant tensors may be carried out in a similar 
manner [15] . For the J\f = 4 theory the energy momentum tensor takes the form 

IV = -V vr+m Eg) + Ai Eg) + A 2 Eg + A 3 Eg] + A 4 Eg) + A 5 Eg) + A a Eg) + A 2 Eg) , (18) 



2 Note that all the terms in (|14p arc of order two in the derivative expansion. Their superscripts are simply 
a means to enumerate them. 
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where 



y(o) 



u x d x a au) + -a au d x u x 

A y(2) _ _ 



si 3 ) 



/' 



y(4) _ ,9, ^ x ^ 



y(i) - /), /i/ . 



si 2 ,) = £/„^a 



(19) 



In the rest of this work, we use the AdS/CFT correspondence to compute the various 
transport coefficients associated with the energy momentum tensor and i?-charge current of 
the M = 4 theoryjfl Sections [3] and @] describe this computation in detail and the results are 
summarized below. 

The energy density and charge density are given by 



3iV 2 

8vr 2 6 4 
/ir 2 iV 2 



where TV is the rank of the gauge group and 



(20a) 
(20b) 



r + = ir 



TlT 

2 



7T 4 T 4 



2 u 2 
1 + -- 



3tt 2 T 2 



2 u 
1 + - 



3tt 2 T 2 



+ 1 



2 u 
1 + ^ 



3tt 2 T 2 



:21a) 



(21b) 



Later we will also need 



1 2 

2 + 



\ 



in 



2/j 2 

37T 2 T 2 



(21c) 



3 In practice our results can be generalized to any CFT whose dual can be truncated to Einstein-Maxwell 
theory on AdSs with a Chern-Simons term. The only difference between the transport coefficients of that 
CFT and the corresponding ones in the AT = 4 theory is an overall multiplicative factor of order unity 
associated with the volume of the compact manifold. 
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The first order transport coefficients are given by 

rj 1 

s An 

- = -r 7 ,Tb 8 
X 2 + 

Q fi 2 r 4 + b 4 
X 



2v / 37r 2 T 2 ' 

where s is the entropy density and x is the susceptibility, 

#2 3 



1 de 
3df 



2vr 



dp 



-N 2 T 2 
A 



The thirteen second order transport coefficients are 



c 

Ai 
c 

As 
c 

A3 
c 

A4 

c 

A5 

c 

Ai 

c 

A3 

c 



1 + J& 



+ 



1+J& 

6r1 



72 7 r 2 r 2 3 4 (2r 2 + r 2 ) 
2 




In 



rl + 2r1 



277T 2 

-1 +31n(2) 

2167T 4 



216 



r + fj, 3 b 4 
54^ 



7T- 



1 + 



6H 



(22a) 
(22b) 

(22c) 



(23) 



(24a) 
(24b) 

(24c) 

(24d) 

(24e) 
(24f) 
(24g) 
(24h) 
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and 



6 _ ln(2) f 1 + o( E U (24i) 



c 72vr 4 T V VT 
^ = (r® + 2rtr\ - / - 2r\% + 4r^) (24j) 

(241) 



c 48vr 2 V 6r 2 
£ 4 /i 3 r+5 9 / /i 2 



c 108tt 2 V 6r 2 
& = /i 2 r^ 8 / _^ 
c " 72 v / 37r 2 V 6r + 



1 + £r h (24m) 



where we defined 



Q2 Q N 2 r 2 

Strictly speaking, all our expressions for the transport coefficients are valid when fi/T 
is not too large. In terms of the gravity dual from which these coefficients were obtained, 
this implies that the black hole is not close to extremality. For brevity, we have omitted the 
somewhat long expressions for A4 and £1 and included only their leading contribution when 
expanded in a small fi/T expansion. The interested reader is referred to appendix [B] for the 
full expressions. A discussion of these results can be found in section [51 We mention here 
that the \x — > limit of t w , Ai, A2 and A3 coincides with the computation of [151 HH EB] f° r 
the n = case. In [16] the dispersion relations for the current were computed via the Kubo 
formula and compared to the Israel- Stewart theory. 



3 Setup 

In the previous section we have explained how one can define the hydrodynamic transport 
coefficients of a conformal theory from the form of the energy momentum tensor and current. 
In what follows we explain how these transport coefficients can be computed from the bulk 
dual of the gauge theory. 

Our starting point is the five dimensional action of Einstein-Maxwell theory 



S 



1 



16ttG, 



g [R+ 12 



4 



1 



12-S/3 



MNOPQ A P P 



d 5 x . 



(26) 
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The metric 



ds 2 = -r 2 f(r)u^u u dx^dx u + r 2 'F ' flu dx )J 'dx u - lu^dx^dr 



(27) 



with m m a fixed four vector satisfying u^u^ 



1- 



f(r) = 1 + 



Q 2 



1 



(28) 
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and P^v as in (jUJ), together with the gauge field 



A r — , 



(29) 



are solutions to the Einstein-Maxwell equations 




Rain + ^9mn 




(30) 



derived from (126]) rl In what follows, Greek indices run over the boundary coordinates [i = 
0, . . . , 3 while Roman indices run over the bulk coordinates N = 0, . . . , 4. 

The metric (12 7p is nothing but a boosted version of the charged black brane solution 
expressed in the Eddington-Finkelstein coordinate system. We would like to extend the 
solution (1271) and (1291) by allowing Q and b to vary slowly with the space-time coordinates. 
It is a simple exercise to check that (1271) and (1291) are no longer solutions to the Einstein- 
Maxwell equations (13"U|) once the boost parameters, charge, and mass of the black hole are 
allowed to vary. Thus, we need to correct the metric (T2T1) and (T2"9"|) to take into account the 
change in w M , b and Q. We will do this order by order in a derivative expansion. To set the 
stage for our perturbative expansion we decompose our metric and gauge field into scalars, 
vectors and tensors with respect to the local velocity field of the fluid, 



ds 2 = r 2 k{r)u p u v dx il, dx v + ^h^P^dx^dx" + r 2 TT iJiU (r)dx^dx u 

+ r 2 j a (r) (P°u u + P> M ) dx^dx u - 2S(r)u^dx' I dr 

= ^g^dx^dx" - 2S(r)u^dx> 1 dr , (31) 



4 Note that we do not require the gauge field to vanish at the future horizon and therefore it is likely that 
it diverges at the past horizon. In fact, it is likely that the whole perturbative solution diverges at the past 
horizon because generic solutions of viscous fluid dynamics are not expected to be regular in the infinite 
past. We thank A. Karch, D. Son, and especially R. Loganayagam for clarifying this point. 
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and 

A r = , A^ = a u {r)P; + c(r)u^ . (32) 

We point out that the various functions k(r), h(r), etc. are not only functions of the radial 
coordinate r but could also depend, in principle, on the charge, Q(x a ), the mass parameter 
b(x a ) and the velocity field or on their derivatives. So k(r), h(r), ... implicitly 

depend on the transverse coordinates. We have chosen an axial gauge for the gauge field, 
A r = 0, and set g^ r oc with g rr = 0. There is one extra gauge degree of freedom which 
we will fix shortly. Our goal is to compute the functions k, h, ir^, j a , S, a v and c order 
by order in a derivative expansion of Q(x a ), b(x a ) and m''(x q ). Using a superscript (n) to 
denote the n'th order contribution to such an expansion, we can rewrite (1271) and 029|) as 

fc(°)(r) = -f(r) , S (0) (r) = l, h<®(r) = 1 , 

j(°)(r) = 0, vr(°J(r) = 0, (33) 

c (0) W = -^, aW(r) = 0. 

A computation of the various functions k(r), h(r), ... is carried out in section HI Once 
these are obtained, we can compute the energy momentum tensor (T^) and i?-charged 
current (J M ) of the boundary theory using the standard AdS/CFT dictionary [3H |35| 1361 [37] 
[381 [39] adopted to Eddington-Finkelstein coordinates [18]. In the Landau gauge this reads 

16ttG 5 (T^) =k® {Au^u v + + 4ttJ5 

(J") =^=^S Tcn [Af\g®] = ~^r^af\ (34) 
>/-g®8Af 8?rG 5 

where a barred superscript (n) indicates the n'th term in a large r (near boundary) expansion 
of the appropriate expression. Working in the Landau frame also requires that 

j&n) = o k^ n) = c (5 ' n) = (35) 

for n > 1. If (135j) is not satisfied this would correspond to a small shift in the local velocity 
fields, the energy and charge densities. 

The chemical potential u of the boundary theory is given by the difference between 
the value of the temporal component of the gauge field at the horizon and its value at the 
boundary (of the unboosted black hole solution). The temperature T of the boundary theory 
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can be obtained from the Hawking temperature, 

M = ^ + >-^) = ±f, T^-^)'-^). (36) 

We have defined r + to be the larger of the two positive roots of f(r) and r_ the smaller of 
its two positive roots (the other four roots of /(r) are given by — r+, — r_ and ±zy/r+ + ri). 
By manipulating ([36]) and f[28l) one can obtain ([21]). Note that an extremal black hole is 
obtained when r + = r_. In the boundary theory this corresponds to the limit T — > 0, while 
keeping fi ^ 0, as can be seen from (12 left . In [ID] it was shown that in this limit, even 
though the temperature vanishes, the mean free path is not necessarily vanishing and one 
might expect a hydrodynamic description of the theory in this regime. However, following 
[JT], it was also argued that this regime of the theory is likely to be unstable. We will also 
see shortly that, from a bulk point of view, our perturbative analysis breaks down when the 
black hole is close to extremality. Hence, in what follows we will assume that fi/T <C 1. 

It is now straightforward to derive the energy density e in (I20al) and the charge density 
p in f|20bl) by inserting f[33]) into f[3~4"[) and using 

ift n 8Vo15 an\ 
where V0I5 = 7r 3 for the M = 4 theory. 

4 The derivative expansion 

If u^, Q and b are constants then (I3"3"|) is a solution to (1301) . As emphasized earlier, if we allow 
the fields w M , b and Q to vary with the space-time coordinates then (1331) is no longer a solution 
to the equations of motion. However, if we allow w M , b and Q to vary slowly in the transverse 
coordinates then we can construct a solution perturbatively. At zero order we have the 
solution ( 1331) . At first order, we look for a correction to ( 1331) . expressed in terms of functions 
fcW(r), hW(r), 7r$(r), j£\r), 5«(r), a { J\r) and c 1 - 1 ) (r) which depend on one derivative of 
the hydrodynamic fields. That is, we insert f[3~ll) and ([32]) . with k{r) = fc^(r) + fc^(r), 
/i(r) = h,(°'(r) + /iW(r), etc., into the Einstein-Maxwell equations f[30]l . omitting all terms 
which contain two or more derivatives of the charge, temperature or velocity fields. These 
equations will of course be linear in fcW(r), h^(r), etc. If these equations can be solved then 
the solution will give us the metric and gauge field of a charged AdSs black hole, where the 
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charge, mass and boost parameters slowly vary in the transverse coordinates, valid to first 
order in gradients of these parameters. With the first order solutions at hand, this procedure 
may be repeated to obtain fc^(r), h^ 2 \r), etc. — a solution to the Einstein-Maxwell equations 
involving two derivatives of the hydrodynamic fields. Up to some caveats which we discuss 
below, one may, in principle, carry out this algorithm to an arbitrary order in the derivative 
expansion. 

It is straightforward, but tedious to compute the order n Einstein-Maxwell equations. A 
method developed in [T3] which simplifies this task is to consider the equations of motion 
in the neighborhood of a point Xq but at arbitrary radial coordinate r. The hydrodynamic 
fields are expanded in a Taylor series around Xq up to order n. Thus, no information is lost 
regarding an order n derivative expansion. Once a solution is obtained around Xq it can be 
uniquely extended to the entire manifold. The interested reader is referred to [H] for an 
extended discussion of this method. As in [H] we choose Xq = 0. At this point we can also 
choose u M = (1, 0, 0, 0), b = b and Q = Q . 

After implementing this technique, we find that in the neighborhood of Xq the Einstein- 
Maxwell equations (1301) take the form 



d r (rVW*^) 


= P^ir) 


(38a) 


d r (r%t\r) + 2V3Q a ( ?\r)) 


= *t\r) 


(38b) 


3<9 r S (n) (r) - \r~ l d r (r 2 d r h {n \r)) 


= S (n) (r) 


(38c) 


d r (r 4 fc (n) (r)) + 8r 3 S (n) (r) + b^ 4 (l - 3r%) d r h {n) (r) - 4=Qo<9 r c (n) 

v3 


= K (n) (r), 


(38d) 



and 



d r (r 3 d r c (n) ) - 2V3Q d r S in) + 3V3Q d r h {n) = C (n) (r) (39a) 
d r [r z f{r)draf\r) + 2v / 3Lg jf >(r)) = A, (n) (r) (39b) 

at order n in a derivative expansion. In addition, there are four constraint equations which 
restrict the allowed values of Q(x a ), b(x a ) and M^(x a ) and reduce to the conservation equa- 
tions ([2]) when expanded to order n. 

While the "kinetic" terms for the unknown fields k^(r), h^ n \r), etc are identical for all 
n, the source terms, on the right hand side of fl38l) and (1391) must be determined at every 
order. We obtain the explicit form of the n — 1 and n = 2 sources in the next section. Once 
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the sources are known, it is simply a manner of integrating the equations of motion (|38j) and 
(I35j) to obtain a solution at order n. 
For the tensor modes irffl we find 

where and reduce to 7r,-™^ and of (138aj) when expanded around Xq = 0. The 
boundary conditions we have imposed are that the boundary metric remain flat, i.e., the 
bulk metric is not deformed near the boundary, and that all singularities are veiled behind 
the outer horizon roughly located at r = r + . The upper limit of the outer integral in (|40p 
ensures that the former boundary condition is satisfied. The lower limit of the inner integral 
in ( 1401) ensures that the outer integrand remains finite at r = r + where /(r+) =0. As 
is standard for charged black holes, once the solution ( 133"]) is perturbed, the inner horizon, 
located at r ~ r_, becomes singular (see for example [32]). As long as this singularity is 
located behind the outer horizon, we should not worry about this. However, since the outer 
horizon is no longer located precisely at r = r + , and since we do not want the fluctuations 
of the horizon to reveal the singularity at r = r_, we require that r + ^> r_. More details 
about the geometry of the perturbed horizon (in the uncharged case) can be found in [33] 

In order to determine the energy momentum tensor of the boundary theory, we do not 
need ^^(r) but only its fourth order term in a near boundary expansion tt^J, c.f., f[3~4"|) . 
From (1401) one finds 

4tt (I) = lim 

au 

^ r— >oc 

Thus, if it is only the order n boundary theory energy momentum tensor we are looking for, 
^tu n \ we are excused from doing the double integral in (HP]) . 

Before considering the integral solutions for the scalar modes, k^ n '(r), S^ n \r), c^ n '(r) and 
h^ n '(r), we recall that we have not completely fixed our gauge for the metric. As in [18], we 
choose the gauge h(r) = 1 since it allows us to easily decouple k^ n \r), S^(r) and c^ n \r). 
Other possible gauges are S(r) = 1 and S(r) = —3/2h(r), which were used in [131 [33]. After 



m=0 



-1) 



rn+1 f\m 



OT^(r) 



[m 



1)! 



P lMV (x)dx 



(41) 
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S (n \x)d 



choosing the gauge h{r) = 1, we find 

S {n) (r) = - 
c (n) (r) = - 



X 



X 



C (n) i 



X 



7S QS 



X) 



8r 3 S {n \x) + -^=Q<9 r c (n) (x 



dx'dx + 2 

+ C r~ 



(42) 



where we required again that the boundary metric is flat and that there are no singularities 
for r > r+. The sources S, K and C reduce to the ones in ( 138c[) . ( I38d|) and ( 139a[) when 
expanded around Xq. The extra integration constants Cq and Co are fixed by the Landau 
gauge (1351) . As we discussed in sections [2] and [31 the charge density and energy density were 
defined relative to the rest frame of a fluid element. This definition was a result of our choice 
of frame. In the gravity dual this choice manifests itself as a choice of adding an extra zero 
momentum quasi-normal mode to c^ n \r) and k^ n \r). Choosing the Landau frame implies 
choosing c^ 2,n ^ = and /c( 4>n ) = in the bulk theory for n > 1. Thus, we choose values for 
Co and Cq such that the second term in a near boundary expansion of c^ (r) and the fourth 
order term in a near boundary expansion of k^ n \r) vanish for n > 1. We refer the reader to 
PH QI] for details. 

For the vector equations we find 




■H(r) 



H 



X 



x3 fW I dx 



C v H x {r) + D v H 2 



(43) 



after extending the solution to the entire manifold. The C u 's are chosen so that the r -4 
terms in a near boundary expansion of j(r) v vanish — again a result of our working in the 
Landau frame. The coefficients D v are chosen so that o^\r) and jv(r) will be finite at the 
horizon. The columns of the matrix H(r) are given by the solutions to the homogeneous 
version of the vector equations, Hi(r) and H 2 (r). Its explicit form and a detailed discussion 
of the solution (|43p and its derivation can be found in appendix As was the case for the 
tensor modes, in order to compute the order n contribution to the i?-charge current we do 
not need the full solution a^\r) but only its r 2 coefficient in a near boundary expansion. 
From 0431) we find 



,( 2 ) 



lim 

r— >oo 



- ( rA u (r) 



- s/3b 4 Q C v - ^-b A Q D v 



(44) 
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where C v and D v are given by 



m=0 + 

and 

A, = -V3Q f° -if + Q 2 r • (46) 

J r+ x 2 ¥ J r + x A J r+ x b 

One caveat in the previous analysis is that we have assumed that a solution exists, i.e., 
that all the integrals are well defined. However, if the sources are too divergent near the 
boundary we would find that the outer integrals in (140]) . (I42p and (143|) do not exist. This 
would imply that one could not impose the boundary condition that the boundary metric 
remains flat. By inspection, one can check that the order n sources for the metric and gauge 
field will not deform the boundary theory as long as 

pW(r) = C(r 2 ), S (n) (r) = C(r" 2 ), K (n) (r) = 0(r 2 ) 
jH(r) = C(r 2 ), C (n) (r) = O(r ), A^(r) = O(r ). (47) 

We have checked that this is the case up to order n = 2. 



4CV = - lim 

r— >oo 



4.1 First order expansion 

It remains to evaluate the various sources, P^, S, K, 3 U , C, and A„, order by order in a 
derivative expansion. By direct computation, the first order sources are given by 



P$(r) = "3rV ( 48a ) 

S (1) (r) = (48b) 

K (1) (r) = 2r 2 <9 / X t (48c) 

jW(r) = 3rV9 M M I/ (48d) 

C (1) (r) = (48e) 

A«(0 = ~ (i?0 a Q + Qu"^) + , (48f) 



when restricted to the neighborhood of Xq. After verifying that (j4"Tj) holds, we can use (jUJ) 
to compute the r~ 2 coefficient of a^\r). Before doing so, we make a few remarks on the 
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index structure of the vector modes. Under the conservation law (|2j), the source terms for 
the vectors can be rewritten as 



jW(r) =J D » vPdpUv 

A«(r) =A Du (r) u^u v + A K (r) P^ Q f + A n (r) 4 , (49) 



where 



A^( r ) = j^( r ) = 3r 2 (50) 

A"(r) = ^ A» = -^ T . (51) 

'• 2 (l + ^) 

Since the differential equations (]38p and ([3"9"P are linear differential equations in the radial 
variable r, the index structure of the sources in the transverse dimensions will carry through 
to the gauge field a^\r). Thus, to first order in the derivative expansion, we already see 
from (Pi9"]) that may be decomposed into terms proportional to £ v , P"d a ^ and u a d a u v . 
Since ay is proportional to the boundary current (J M ) via (I34p . and the boundary current 
is Weyl invariant, this implies that the non-Weyl invariant term u a d a u u can not contribute 
to a„ . Indeed, from (jUJ) we find 

4 5) = -V3Qu„ + ±-K 2 r 7 + T*b 8 d^ - ^Q 2 b%. (52) 

Using (1341) we see that the expectation value of the -R-charge current ( J At ) takes the form 
(II 7p with Q and k as in fl22cD and fl22bl) . The energy momentum tensor can be evaluated in a 
similar manner. We compute irffl from (jHJ), and from (j34j) we obtain (T^). Not surprisingly, 
it takes the form ffl~8]) with a shear viscosity 77 as in (122 ap . 

Had we been content with the first derivative corrections to the energy momentum tensor 
and current, we could have stopped here. Since we will be computing also the second order 
corrections, we need the full gravity solution to first order in a derivative expansion. Almost 
all the sources for the scalar terms in the metric are trivial, and integrating them gives us 

S<U(r) = c (1) (r) = fc«(r) = —d,^ . (53) 

3r 



For the tensor modes we find 



*<$(r) = F(r)a ia/ , ( 54 ) 
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where 



F(r) = V Y (r ''" 3) 4 ln(r ~ X) = 1 - 4 + °(r~ 5 ) (55) 



and the sum runs over all six roots of fir) which was defined in ( l28l) . i.e. ±r + ,±r_ and 
±iyV+ + r^. Note that F(r) depends implicitly on the transverse coordinates. The expres- 
sions for the vector components of the metric are somewhat more involved. By explicitly 
carrying out the integral in (fl3l) we find 

tiXr) = - -u^u v - b ^t+r(r)d^ 

a«(r)=^(r)9^-^^4, (56) 

where a K {r) and j K (r) are given by rather long expressions which we avoid writing out 
explicitly here and which are given in appendix [Bj equations ( 1821) and ( |83i) . Their leading 
order behavior is given by 

2 7 T s b 8 

f{r) = 0(r- 6 ) . (57) 



4.2 Second order 

At second order we find that the scalar sources for the Einstein-Maxwell equations (1391) and 
(138]) satisfy (|47|l as required. The sources for the tensor modes are given by 

P<2 = P r 'EW + pAi S W + P * E g) + pAs E (3) + pA 4S (4) + pA 5E (5) + pAxgjl) + pA 2 g(2) (5g) 
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with 



P T '(r) 
P Al (r) 
P A2 (r) 

P A3 (r) 
P A4 frl 



P As (r) 
P Al (r) 

P A2 (r) 



r - (r 3 F(r))' -r 3 F'(r) 

r - 3r 2 F(r) - (r 3 - r 3 )F'(r) 

2 fr - (r 3 F(r))' - r 3 F'(r)) - 4r 
1 2Q 2 \ 4Q 4 6 4 



4r 1 



fo 4 r 4 



3 (ttT) 2 + r 2 - lnTr 
~2 T tx (r 2 + (ttT) 2 ) 



-3fc 4 + 4 - 



12Q 2 6 4 



7T 



T 3 3r 2 \ 



+ 1 4^+4^;' ln 

-2 (r 3 r(r))' 



(ttT)* 



(r + irTf 



2 arctan f ) 



7T 



v^^Wr 4 v^g 2 6 4 f 12Q 2 a 

-a (r) - 2r f(r)a (r) 



1 + ^ ) r 4 



6r; 



6Q 3 6 4 



0*) 



(59a) 
(59b) 
(59c) 

(59d) 



(59e) 
(59f) 
(59g) 

(59h) 



The coefficients A, and tv of the energy momentum tensor can be computed by inserting the 
corresponding sources of (1591) into (14~T!) and using (1341) . The results are listed in (1241) . The 
full expressions for P A4 (r) can be found in (I80aj) . 
The sources for the vector modes are given by 



A (2) _Aftw(l) , a6t?(2) , A&w(3) , a£i^(4) , a|itj(1) 



(60) 



where 



A 6 (r 
A 6 (r 

A 6 (r 
A ?4 (r 
AWr 



ra K, (r) - T%*ii 2 r 7 + F' (r) 
4r +/ u 2 



7r 2 r 7 T 3 6 8 



1 



3r3fl + ^ 







4x/36 5 Q 3 2v/36Q 



V3Q S 



F(r)V b A 
+ 2H 



6 4 F'fr N 



(61a) 
(61b) 

(61c) 
(61d) 
(61e) 
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and 



J ?1 (r) 




(3r 5 + 4vrTr 4 + 7vr 2 T 2 r 3 + 10vr 3 T 3 r 2 + 6vr 4 T 4 r + 4vr 5 T 5 ) //r 



+ O (fi 2 ) (62a) 




J ?3 (r) 




(62b) 



(62c) 



J« 4 (r) 



26r 



(62d) 





) 



(62e) 



The full expression for Jfr(r) can be found in (J80b|) . The coefficients & and & in (123501 
were computed from (l6Tj) and (1621) with the help of (1441) and (134|) . 

5 Discussion 

To summarize, on the gravity side our results show that one can extend the Reissner- 
Nordstrom AdSs black hole solutions to black holes for which the charge density varies with 
the space-time coordinates. The field theory dual of this configuration is a conformal fluid 
with a non vanishing chemical potential and non trivial i?-charge current. The transport 
coefficients of this fluid are listed in (J20l) . (J22l) . (J24l) and ( |79l) . Once the chemical potential 
vanishes one obtains the transport coefficients calculated in [T4"l [T5l [16] for the CFT dual of 
an uncharged black hole. 

At first order in the derivative expansion we have found a transport coefficient of the 
i?-charge current, Q, associated with the vorticity (see (jHJ) and (I22c|l ) which, as far as 
we know, has not appeared in the literature so far. One can trace back the appearance of 
this component to the Chern-Simons term in (1261) . If the Chern-Simons term were absent 
from the Lagrangian, Q would vanish. A similar statement can be made for the second order 
transport coefficients fi, Ai and A 2 E 

Our result for the shear viscosity to entropy ratio agrees with ([T]). This observation has 
already been made in [5], [231 121] and seems to be a feature of any gauge theory with a 
holographic dual. In light of the universality of r)/s, it is natural to inquire if there are other 

J Notice, however, that the presence of the Chern-Simons term is required by supcrsymmetry. 
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hydrodynamic quantities whose value is universal [131 HH 05] . In [15] it was suggested that 
the ratio of the electrical conductivity to the susceptibility of certain CFT's might have a 
universal value in the same sense as (JT]). In the \i — > limit our result for k/x m f l22bl) 
is consistent with the prediction of [15]. In that context it is also interesting to compare 
our full expression for k/x in f!22bj) to the one obtained from the results of [5], where the 
chemical potential of a different U(l) subgroup of the 50(6) i?-symmetry group was turned 
on. As far as we can tell, the ratios are rather different. Also, the constant of proportionality 
in the analogue of the Wiedemann- Franz law [46J, in our case 



differs from the one obtained in [5] by a factor of two. 

We have normalized all the second order transport coefficients (T51]) relative to the second 
derivative of the energy density with respect to the temperature. This was done in order to 
conform to [IB] and in order to get rid of some multiplicative constants related to the number 
of degrees of freedom of the gauge theory. Clearly, the expressions in (1241) are different from 
their fi = counterparts. However, it is interesting to note that 



for any value of and that this relation also holds, in the n = case, in more than d=4 
transverse dimensions [T8] . It would be nice to check if this relation remains valid even when 
fi 7^ and d > 4, or in other theories with a holographic dual. 

Another interesting transport coefficient is A3 which vanishes when /i = both at strong 
and weak coupling [15]. When fi 7^ it does not vanish, at least not at strong coupling. 
From a gravitational point of view, one reason for this difference is the Chern-Simons term 
in fl2"B"|) . As discussed above, this term is responsible for the contributions proportional to 
in jji (r) and a^\r). These expressions contribute to the source P^(r) via terms quadratic 
in jji (r) and a^\r). One can verify that this is the origin of the second summand on the 
right hand side of (159dl) . However, even in the absence of the Chern-Simons term A3 would 
not vanish when /1 7^ 0. Of course, the choice of basis in (ir§j) is not uniquely determined, it 
is always possible to redefine 




(63) 



4Ai + A 2 = 2r/r^ 



(64) 




+ 



(65) 
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so that the coefficient of S 3 vanishes both in the \x 7^ case and in the \x — > limit. Further 
physical guidance would be needed in order to decide whether a certain choice of basis is 
preferred over another. 



Appendix 

A The vector modes' equations of motion 

In section H] we have shown that the equations of motion for the perturbations of the metric 
and gauge field (!30|) involve two vector modes ffl{r) and a^\r) which are coupled, 

d r (r 3 f{r)d r a{r) + 2V3Q j{r)\ = A(r) 

d r (r 5 d r j{r) + 2V3Q a{r)) = J(r). (66) 

We have removed the superscript (n) specifying the order of the solution and the vector 
subscript fi for clarity. In this section we show how to obtain the solution to these equations 
for arbitrary sources A(r) and J(r). The result was already anticipated in (143|) . 

We start by considering the homogeneous version of these equations. The four homoge- 
neous solutions are given by 



<r)\ = l\ a(r)\ = 
J(r) Ur \j(r) \l 




V3Q b* 

1 _ MQo 

r-4 rpt 



a(r)\ = / -^ + ^i^a(x)ln|r 2 -x 2 | 

j(r) J \\b* + ^ + 2V3Q f(r) Z x a(x) In |r 2 

where we have introduced the following notation: 



x 2 \ 



= #i(r), (67a) 
# 2 (r) , (67b) 



4(-4 + 276^)' 2 

-2^P Qo/M , , v^, 2 Qo(x 2 + 3^Q 2 



61 j P 2 (r) ' 8(-l + 3x^)(-4 + 27TO) 

and x runs over the roots of /(r) which we denote by ±r + , ±r_ and iz^/r^ + r 2 with 
r+ > r_ > 0. The outer horizon of the unperturbed black hole is located at r = r + . The 
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two constant modes correspond to a deformation of the boundary, Hi is a zero momentum 
quasi-normal mode corresponding to a shift in the boost parameters of the black brane, and 
H 2 is a solution which diverges at the horizons. In particular, at the outer horizon we find 

lim H 2 (r) = 2(2^F 1 +' ^ +M ( 69 ) 

The rest of the solutions are finite at r = r+. 

To solve the non-homogeneous equations we integrate (I66I) once, 

«9 r a(r) + 2v / 3Lg 4^T = ^WT ( 70a ) 

/r 
3(x)dx. (70b) 

For the moment we keep the lower limits of integration unspecified. The solutions to the 
homogeneous version of flTOj) which are first order equations, can be obtained from linear 
combinations of (1671) . They are given by the columns of 

( ^P- -^ + ^#V r a(x)ln|r 2 -2 2 | 

H(r) =\ ' P ( \ ,— . 71 

\2V3Q f(r) ^ + 2V3Q f(r) J2 X a(x) In |r 2 - x°' 



The overall multiplicative factor of the homogeneous solutions has been chosen so that 
\H\ = I. 

With the homogeneous solutions at hand, one can use the method of variation of param- 
eters to find a particular solution to the inhomogeneous first order equations (jTDj) . This is 
given by 



a(r) 



, pA{x')dx' . 

-H{r) H-\x)[ _*: f y^,)dx. (72) 
j(r) J J r \x 5 J 3(x')dx J 

We choose the integration constants by requiring that the metric is differentiable up to 
and including the outer horizon and that there is no deformation of the boundary metric. 
Since H is finite at the boundary, the latter requirement implies that we should set the upper 
limit of the outer integral to infinity. We assume this integral exists, i.e., 

Jto(r) = C(r 2 ), AW(r) = O(r ). (73) 

If (1731) does not hold then there is no asymptotically AdS solution. As stated in the main 
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text, we have checked that (!73|) is satisfied up to second order in the derivative expansion. 
The other requirement, that the metric is differentiable, implies that we should set the lower 
limit of the inner integral to r + and add an appropriate multiple of the homogeneous solution 
H 2- Indeed, once the lower limit of the inner integral is set to r = r + then the outer integrand 
in (1721) will be finite at the horizon. Thus, the only terms which may diverge at the horizon 
can arise from the logarithmic divergence in H multiplying the outer integral in (1721) . Since 




-H(r) 



f* A(x')dx' 



VZQo 

2 \2 



2(27rTr 



H 



In \r 



f x J{x')dx' 




dx 



6Ql 



x- 



J 3(x')dx'j dx 
+ O(r + ), (74) 



then according to (1691) one can get rid of the remaining logarithmic divergence by adding a 
term proportional to H2(r) to our solution. We are still left with one integration constant: 
the homogeneous solution Hi neither deforms the boundary nor diverges at the horizon, 
so we may add it to (172]) without spoiling the boundary conditions. As we have mentioned 
earlier, Hi is the homogeneous solution associated with a shift in the boost parameters. From 
the point of view of the boundary theory, this corresponds to an ambiguity in the definition 
of the velocity field which is fixed by going to the Landau frame. Fixing the Landau gauge in 
the boundary theory corresponds to setting the fourth order coefficient of a near boundary 
expansion of j( n )(r) to zero, cf. (1351) . This precisely fixes the remaining integration constant. 
Our final result for the solution to (1661) is then 




j; + A(x')dx' 

H'Hx) I x3 fW I dx + CHiir) + DH 2 (r) 
i" 5 f* J(x')dx' j 



(75) 



C = - 



lim 

r— >oo 



m r m-\-l 




d™J(r) 



J(x)dx 



vm=0,l 



32m+l 



A(x)dx 



(76) 
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by requiring that = and 

by requiring that the metric is differentiable at the outer horizon (notice that this does not 
interfere with demanding = since H2 does not contain any term proportional to r -4 ). 
We point out that since A = O(r ), c.f. f J73|) . the second line in (1761) will always evaluate 
to zero. Also, in practice it is efficient to replace the double integrals in ( 1771) with single 
integrals. This can be done by integrating by parts. Using (1751) . we obtain 

D = -V3Q ^-dx -if ^rdx + Ql ^-dx . (78) 

Jr+ % % Jr+ % Jr+ % 

Equation (|4"4"1) in the main text was obtained by expanding ( 1751) in a series expansion 
near the boundary, using (1751) and extending the solution from the neighborhood of Xq to 
R 3 ' 1 . 



B Long expressions 

In certain places in the main text the expressions we have found were somewhat long. In 

this appendix we have collected the expressions which were omitted. 

The full expressions for the transport coefficients A4 and ^ whose expansion (in ^) 
appeared in section [J2, equations (I24el) and (I24il) . are given by 

A 4 (Sr 1 ! 3 + 36r\r]} + 372r 4 _rL 2 + 990r5j_ri° + 1523r**_ r » + 1438r^°r1 + 696r| 2 r* + 136r]^r 2 _ - 16r}_ 6 ) T 4 
3456r 4 (\+ (2r 2 _ + r 2 ) 2 (r 4 + r\ r 2 _ + r± ) 3 

(^+"^-^) 6 ( 1 +4) 1 "(^fe|) 



(79a) 



and 



8tt 2 r+ (r 2 + r 2 ) (l3ri + 13r 2 r 2 + 10r 4 ) (ri + r 2 r 3 - 2r%rJ)' 
Jp^ 1 ~ 32vr (ri + r 2 ri + r %f (ri + r 2 ri + 2r 4 ) 



(ri + r 2 ri — 2r 4 ) 3 . / ri — r 2 



!>7rr + (ri + r 2 ri + r 4 ) (2ri + 3r+ri + 5r 4 r 2 + 2r+) \2r+ + ?" 



^ ttV-t • (79b) 
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The source terms missing from (I59el) and (I62al) in section [4] are 



p ^ 4 = y^T^(,i +r ^-2^) 3 

12V3Tvr_r 2 _ (ri° + 3r^r* + 6r^r^ + 7r^r* + 5r* ri + 27-1°) ' 



12^3n-r_r^ (ri° + 3r 2 ri + 6r* r* + 7r^ r 4 + 5r^ri + 2r^_°) 



\3 



5 ( r ~ +r+r " 2r+ ^ r ( 2L^l\ +2 V3r_r + Jrl+rl (a" (r)j K (r)) ' ' + ' " " +; ^ ."'(r) 

127rr + (rS +2r2 r 6 +4rlrt +3r6rl +2rS) V r / V 8 7rr 2 (ri + r 2 ri + r* ) " 

V^( r 4 +r 2 r 2 _ 2r 4)3 ( r 4 +r 2 r 2 _ 2r 4)3 
4- i i 3V r 3 3 ( Q1 W 7 - — r- — (80a) 



(ri + 2r^.ri + 4r* ri + 3r^ri + 2r^ ) ^ ' 4-/rr + (ri + 2r^ ri + 4r* ri + 3r^ ri + 2r^ ) r 2 



and 



Jft(r)= ((r 5 F'(r)r(r))' 



? + + r 2 r_(r - r + ) 2 r^ (3r + (r 2 + 2r + r + 3r 2 ) r 2 + (2r + r+) (ri + r\ r 2 + 4r+)) T 2 

+ 2r 2 r & _ + 4r^ ri + 3r-6 r 2 + 2r% ) f(r) 2 

(80b) 

where we have defined 

d,f{r) = j« d ®(r)u a d aUfl + f [dQ \r)d^Q . (81) 

Finally, in (I56I) we have parameterized the first order bulk solutions for the vector modes 
by two functions a K (r) and j K (r). The first, a K (r), is given by 

jrl + r 2 (ari + 3r 2 ri + 2r j ) (r 2 - ri ) 3 3ri r + (ri + 2r= ri - r* r 2 _ - 2r% ) 

a (r) = 



4r + (2r 2 + r 2 ^) (ri + r 2 r 2 + 2r* ) 2irr(2r 2 + r 2 ) (r£ + r 2 r 2 + 2r* ) 

27r^_ (ri + ) 2 (ri + r^r 2 - 2r* ) ri 3r + (ri - r\ ) 3 ^/ri + r 2 (2ri + 5r+ri + r>r\r 2 _ + 2r^_) ri 



+ 



2wr 



16r 2 7T (2ri + Sr^ri + 3r* ri + r 1 ^) Sr 2 (2ri + r 2 ) (ri + 2r 2 ri + 4r* ri + 3r^ri + 2r^ ) 

(r 2 _ -r 2 ) 3 (2ri +5r 2 ri +5rjri + 2r% ) 3rir + (ri -r 2 ) 3 (ri + r 2 ) 3/2 (2ri + 3 r 2 ri +2rj) \ 

+ yjr 2 _ + r\ {2r 2 _ + r 2 ) 3 (ri + r 2 ri + 2r*) 4,rr 2 (2r£. + r 2 ) 3 (ri +2r 2 ri +4r* ri + 3rS ri + 2r» ) j ^ ^ jr 2 _ + r 2 
(ri + r + r_ + r^) (r'i - r + ri + 2r^_r_ - 2r5|_) 3 3ri (ri + rjj.) (ri - r + ri + 2r^.r_ - 2ri^) 3 



4-n-r 2 ^ (2ri + r^) (ri + r^ ri + 2r* ) 8irr 2 (2ri + r 2 ^)' (ri - r + ri + 2r^ri - r^ri + 3r* ri - 2r^r_ + 2r^) 

r_ (ri - r + r_ + r+ ) (ri + r + ri + 2r^_r_ + 2r^) 3 3ri (ri + r^.) (ri + r + ri + 2r 2 _r _ + 2r^) 3 



ln(r- 



47rr 2 _ (2ri + r^) (ri + r 2 , ri + 2r| ) 87rr 2 (2ri + r 2 ^ ) (ri + r + ri + 2r 2 ri + r^ ri + 3r* ri + 2r^ r _ + 2r^ ) 

( ^f +{r - +Tl) , ~ ~) ^ + ^ + ) + ( ^ll^nL^l + (l^lVl^) ) ln p + r 2 + r 2 
\4^r 2 (ri +r 2 ri + r^) ^ J ^ 8wr 2 (2ri + r= ) (ri + r= ri + ) 4^r 2 (2ri+r 2 ) ^ V 



ln( r 



(82) 
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and j K (r) can be determined from a K (r) through the equation of motion, 

(rt + r 2 r 2 — 2r+) 3 (r (rt + r\r 2 _ + 2r+) — 2r + (rt + r 2 + r 2 _ + r+)) 



16v / 37rrr_r^y / r^ _ +"r^ (ri + r 2 r 2 + r+) 2 (ri + r 2 r 2 + 2r 



_ (r 2 - r 2 ) (r| - r 2 ) (r 2 + r 2 + rjj ^ 
2v^3r 3 r_r + ^/r 2 + r 2 
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